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ABSTRACT 

I develop a detailed empirical model for the chromosphere and wind of 31 Cyg 
based on a previously published analysis of WE spectra from the 1993 eclipse and 
on the thermodynamics of how the wind must be driven. I then use this model to 
interpret observations of single supergiant stars and to assess the evidence that 
their winds are fundamentally different from those of supergiants in the binary 
systems. This model naturally predicts a certain level of clumping of the gas 
to balance the pressure that drives the wind. It also predicts that anisotropic 
turbulence, such as would result from transverse displacements of Alfven waves 
directed along radial magnetic flux lines, would not give the roughly Gaussian 
profiles of emission lines seen in cool giant stars. Furthermore, it implies that 
C II] may not tell us much at all about general conditions in chromospheres. 
Finally, I speculate that chaotic magnetic fields, in dynamical equilibrium with 
the gas of the wind, are the actual driving mechanism. 

Subject headings: stars: winds — stars: chromospheres 



1. INTRODUCTION 

It's no great stretch to assume that all cool stars (later than mid-F spectral type) have 
chromospheres and, likely, winds as well. These are regions above the photosphere too hot 
to be heated radiatively; rather, they must be heated by magnetic processes or by shocks. 
Most analyses such stars' chromospheres and winds are based on interpreting the strengths 
and shapes of emission lines in their spectra. The Sun is the only one for which we can 
actually map out the chromosphere in any appreciable spatial detail. All the others are 
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unresolved, which limits us to interpreting global properties of the emitting gas in terms 
of models without any real spatial information (zero-dimensional). There are a number of 
standard diagnostics used in this work. For cool giant and supergiant stars, we have the 
C II] intersystem multiplet at 2325 A (Stencel et al. 1981), the Al II] multiplet at 2669 A 
(e.g., Judge 1986a, 1986b; Harper 1992), the Mg II h and k lines (e.g., Stencel et al. 1980; 
Milkey et al. 1975), an abundance of Fe II lines in the UV (ultraviolet) affected by scattering 
in winds (Judge & Jordan 1991), Ha (Johnson & Klinglesmith 1965; Mallik 1993), and the 
H I Lyman lines (Landsman & Simon 1993; Haisch et al. 1977). Many of lines involved fall 
in the UV because of the temperatures and densities obtaining in these chromospheres, so 
our knowledge of cool giants' chromospheres necessarily has developed through observations 
from IUE, the International Ultraviolet Explorer satellite (e.g., Judge 1988) and more recently 
HST, the Hubble Space Telescope (e.g., Carpenter et al. cited below). Other observations of 
these chromospheres involve measuring radio thermal emission from their winds (Drake & 
Linsky 1986; Harper et al. 2001). 

An alternative, complementary approach to mapping the properties of chromospheres 
solely by interpreting observations of the emission lines of single stars is to use the ( Aurigae 
binaries (Wright 1959, 1970; Wilson & Abt 1954) to probe their chromospheres and winds 
in greater spatial detail (legitimately at least one-dimensionally). These systems consist of 
a cool supergiant paired with a B dwarf, the prototype system ( Aur consisting of K4 lb 
and B5 V components. Absorptions in the spectrum of the B star by the gas in the wind 
and chromosphere of the K star give us a spatial probe of the K star's wind. Again, the 
telling absorptions are in the UV, so most of our knowledge of these stars has come from IUE 
(e.g., Eaton 1996) with incremental improvements from HST (e.g., Baade et al. 1996). More 
recently Harper et al. (2005) have used radio observations of ( Aur to verify the atmospheric 
structure they derive from UV lines. 

The purpose of this paper is to develop a detailed physical model for the chromo- 
sphere/wind of 31 Cyg, to test it, and then use it to interpret observations of emission 
lines of single stars. There is a commonly held belief, even prejudice, that the chromo- 
spheres/winds of single stars are different than the those of binary components in that winds 
of single stars accelerate much more rapidly than winds of the binaries and reach lower 
terminal velocities. This paper assesses whether the observations of single stars might be 
consistent with the slower acceleration of binary components. I shall begin with a review 
of the picture developed for single stars through traditional analyses (§ 2), then discuss the 
model for cool supergiants developed from observations of the ( Aur binary 31 Cyg (§ 3), 
testing it through comparisons with emission lines and radio continuum emission. I then use 
this model to interpret the measurements for single stars (§ 4) and give some speculations 
about wind acceleration (§ 5). My approach follows the radically different viewpoint that 
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pressures supporting and accelerating the outer atmosphere are a microscopic phenomenon, 
operating on scales smaller than we resolve with our observations, and are not simply the 
manifestation of global waves. 

2. THE GENERAL PICTURE FOR SINGLE STARS 

Our observational understanding of single stars consists of several types of data. (1) 
relative and absolute fluxes of a wide variety of UV and optical chromospheric emission lines, 

(2) profiles of optically-thin emission lines, reflecting the kinematics of gas forming them, 

(3) velocity structure of optically thick emission lines, reflecting bulk flows in the chromo- 
spheres/winds, (4) line strengths of density-sensitive multiplets, namely C II] UV0.01 A2325; 
and (5) strengths of fluorescent lines, reflecting mean intensities in some of the deeper parts 
of the chromosphere. Many of these results come from HST spectra discussed in observa- 
tional papers by K.G. Carpenter and his collaborators (Carpenter et al. 1991,1995,1997,1999; 
Robinson et al 1998); the more imaginative interpretations of them seem to have come mostly 
from P.G. Judge (e.g., Judge 1994). 

The important results of these papers are (1) measurements of electron density in some 
part of the emitting region of the chromosphere, (2) turbulent velocities from optically-thin 
emission lines, (3) possible evidence for multiple components in the profiles of Cll], but 
not other optically-thin lines, as well as in the strengths of fluorescent molecular lines, (4) 
some indication of the temperatures in the chromospheres/winds from excitation of Fe II 
and relative strengths of various emission lines, (5) detection of slight amounts of more 
highly ionized, possibly hotter, gas (e.g., C IV emission), and (6) evidence for atmospheric 
expansion and a crude measurement of the density-velocity profile from self-reversed Fe II 
lines and Mg II h and k. 

Chromospheric Velocity Structure: Many of the weaker intrinsic lines of common 
elements are weak enough in chromospheric spectra as to be optically thin. They should be 
formed in essentially nebular conditions, with every excitation giving a potentially observable 
photon, so that their line profiles would reflect the kinematics and ionization structure of 
the emitting gas. The superior resolution and signal to noise of HST have made it possible 
to measure their profiles reliably. The turbulence measured from such an optically-thin 
line can be conveniently characterized by a velocity, t> =FWHM/1.67, the parameter in the 
Gaussian part of the velocity distribution [exp - ^/" ) 2 )], which Judge (1986a, 1986b) called b. 
Carpenter et al. fit Gaussians to such lines in a variety of cool stars. We may estimate the 
turbulence by averaging their results for roughly nine lines between 1900 and 2850 A, with 
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the C II] multiplet excluded, which give FWHM = 18.2 km s" 1 (a Tau), 29.9 (A Vel), 27.2 
(a Ori), and 23.6 (7 Cru). I have corrected these values for the 13-15 km s _1 resolution of 
the GHRS spectrograph which Carpenter et al. mention occasionally in passing, except in the 
cases they explicitly stated that they had corrected their published values. These FWHM's 
correspond to turbulent velocities in the range vq = 10.9-17.9 km s _1 for the emitting gas. 
This is well supersonic for the temperatures expected in these chromospheres (< 10,000 K). 

The most perplexing result of such kinematic analyses is the existence of broad wings 
and a ~4 km s" 1 redshift of the C II] lines in most of these stars. These may well not be 
formed under similar conditions as most other weak lines, and, in fact, those other lines 
(such as fluorescent Fe I and weak lines of Fe II and Al II]) do not show these phenomena. 
This discrepancy implies some sort of anisotropic turbulence seen only in C II] or, perhaps, 
multiple components in the chromosphere, again not seen in the weak non-C II] lines. Single 
Gaussians fit to the C II] profiles in the various stars give values of FWHM = 24 km s _1 
(a Tau), 28 (7 Dra), 36 (A Vel), 35 (a Ori), and 30 (7 Cru), not corrected for the point 
spread function of the spectrograph. For 7 Cru a two-component fit gives FWHM=27/42 
for core/wings; for a Ori, FWHM=19/48. For A Vel, the star most like 31 Cyg, we have 
a characteristic velocity of 21.6 km s" 1 , presumably uncorrected for the resolution of the 
spectrograph, or 19.9 km s _1 corrected; again, supersonic for the expected chromospheric 
temperatures. Of particular significance would be any differences of width in these optically- 
thin lines; if the C II] lines were formed primarily in a different region from the others, 
they would not necessarily have the same broadening. Indeed, the C II] lines are marginally 
broader than the other lines in a Tau (AFWHM = 1.6 km s _1 ), A Vel (3.4 km s" 1 ), and 
7 Cru (2.7 km s _1 ). These differences must at least partially reflect the non-Gaussian profiles 
of C II]. 

Electron Densities: These come primarily from the density-sensitive line ratios of the 
C II] (Stencel et al. 1981). Electron densities derived from the C II] line ratios in HST and 
IV E spectra of typical stars are near 10 9 cm~ 3 . It is important to remember that these values 
apply only to the parts of the chromosphere/wind where the photons are formed. Because 
C is theoretically expected to be mostly neutral in the deeper reaches of the chromosphere, 
one might expect that the measured electron densities would apply only to the warmer 
parts of the chromosphere /wind. The C II] lines also give limits on the optical depths of 
chromospheric gas. 

Judge (1994) analyzed the variations of very good observations of a Tau with Doppler 
shift to derive electron density as a function of velocity. Again, his results were perplexing. 
They gave higher electron densities, by a factor of four, for positive velocities (away from 
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us; ostensibly toward the star) than for negative velocities. Furthermore, the emission in 
a Tau was shifted to the red by 4 km s _1 , which Judge interprets as a downflow in the denser 
regions of the chromosphere. A redshift of 2-4 km s _1 seems to be ubiquitous in the cool 
giants and supergiants (Judge & Carpenter 1998). On reflection, that the measured electron 
density is so similar in all these stars, in spite of variations in mass-loss rate of several orders 
of magnitude (see Harper 1996), is probably just as strange as the line shapes. 

Line Strengths: Relative strengths of optically-thin lines give a measure of the mass, 
temperature, and electron density of emitting gas averaged over the whole chromosphere. 
They may be interpreted with appropriate empirical or semi-empirical models as Judge 
(1986a, 1986b) did for several bright stars. This sort of analysis provides global constraints 
on the integral 



Fluxes for characteristic lines come from IUE spectra and seem to be roughly consistent 
from star to star among the cool giants. Such collisionally excited lines as Si II] UV0.01, 
Al II] UV1, and Mg II h and k have consistent ratios for cool giants (Judge & Jordan 1991). 
C II] UV0.01 is not necessarily proportional to the other collisionally excited lines (Judge et 
al. 1992); however, it is roughly equal in strength to Al II] for most of these stars. For the 
values given by Judge & Jordan, we can summarize these line strengths as in Table 1. For 
the windy (non-coronal) giants, the flux in C II] A2325.4 seems to be roughly equal to flux 
in Al II] A2669. We have also listed strengths for three ( Aur binaries as best we can judge 
them. The Al II] fluxes for these binaries (Eaton 1992) are relatively well determined. The 
C II] fluxes, on the other hand, are very poorly measured and in a noisy part of the IUE 
spectrum. There are no spectra for 31 Cyg exposed long enough to detect this line, and 
it is hard to differentiate it from the noise in the spectra for ( Aur and 32 Cyg in which 
Schroder et al. (1988) claimed to detect it, although Harper et al. (2005) may have detected 
the strongest component with HST. I have given what I consider the best measurements 
possible for these spectra, but in the following analysis, I shall assume f(A2325.4) « f(A2669) 
in all the windy giants. The fluorescent Fe I UV44 lines can be detected in all three classical 
( Aur systems (e.g., Bauer & Stencel 1989), but these two lines, especially A2844, are highly 
blended with other features. 

Judge (1986a, 1986b; Judge & Jordan 1991) has made simple empirical models for 
three cool giants, finding the emission is probably excited in a gas near 7000 K and an 
electron density of 10 9 (an assumed global value of n e from C II] UV0.01). Carpenter et al. 
(1999) estimated that the Fe II lines observed in the wind of A Vel are scattered by gas at 
~ 6000 K from the relative strengths of various multiplets. Other estimates of the electron 
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temperatures of single stars come from adjusting semi-empirical models (e.g., Kelch et al. 
1978) to give observed line profiles for charactistic emission lines and Ha absorption. 

Fluorescent Lines: The classical fluorescent lines of K giants are Fe I UV44 excited by 
the Mg II k line. These lines are optically thick and must be formed in the deeper, denser 
parts of the chromosphere for there to be any neutral iron to scatter them (Harper 1990). 
Other fluorescent lines are scattered by molecular species (McMurray et al. 1999; McMurray 
& Jordan 2000). These molecular species must be present in some part of the chromosphere, 
but the eclipses of ( Aur binaries generally do not detect them, probably because there is 
plenty of cold dense gas in the inner parts of the chromosphere to hide their absorptions in 
the competing atomic features. McMurray et al. have argued that the strengths of fluorescent 
H 2 and CO lines, which they could not reproduce with their non-dimensional semi-empirical 
model of a Tau, require a multi-component atmosphere, possibly with shocks deep in the 
atmosphere to generate enough Lya flux there to excite these molecules radiatively. 

More Highly Ionized Species: Long exposures in the UV have detected lines of highly 
ionized species even in the windy giants. These include intersystem lines of Si III and C III 
(e.g., Carpenter et al. 1999) and the resonance doublet of C IV (e.g., Robinson et al. 1998). 
In the Sun, such highly ionized species emit lines in the transition region at ~ 50,000 K. 
In the windy giants, this interpretation is problematic because there is no evidence for the 
coronae that create the transition regions in the dwarfs. 

Expansion- Velocity Structure of the Wind: This is very difficult to get from obser- 
vations of single stars. There is ample evidence of winds in the P-Cyg profiles of Mg II and 
from the asymmetric, variable profiles of Ha in cool supergiants (Zarro & Rodgers 1983; 
Mallik 1993). Harper (1996) summarized analyses of wind profiles derived from single stars; 
Harper et al. (2005) have discussed this in somewhat more detail. Much of the evidence 
comes from P-Cyg profiles in the UV. Fe II lines of increasing intinsic strength show increas- 
ingly negative velocities of their self reversals, attributable to wind acceleration, in a number 
of cool supergiants, notably A Vel (Carpenter et al. 1999), 7 Cm (Carpenter et al. (1995), 
a Ori (Carpenter & Robinson 1997), and both a Tau and 7 Dra (Robinson et al. (1998). Of 
these, A Vel had the most extensive coverage of the accelerating wind, with centers of the 
shell lines shifted to ~ —32 km s _1 with respect to the star, and Carpenter and Robinson 
have interpreted the data for these stars to infer mass-loss rates, terminal velocities, level of 
turbulence, and density structure. Such analyses give winds accelerating much faster than 
those of the ( Aur binaries or, for that matter, of single stars as deduced from radio observa- 
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tions (Harper et al. 2001, 2005; Carpenter et al. 1999). In fact, the analysis of thermal radio 
emission has become a fruitful technique for deducing the wind structure of single windy 
giants (e.g., Drake & Linsky 1986; Harper et al. 2001). 

3. THE MODEL FOR 31 CYG 

Although optical spectra provided many insights into the nature of the extended at- 
mospheres of the cool supergiant components of ( Aur systems, the field was essentially 
dormant from the mid 1950's until IUE revitalized it with panchromatic UV spectra. I have 
discussed the results of such studies in a review at the Cool Stars 9 meeting (Eaton 1996). I 
shall summarize them here as follows: First, the new UV observations recorded absorptions 
from most of the important species expected to exist in chromospheres and winds of cool 
giant stars. In contrast, optical spectra give very few of these species. Furthermore these 
UV absorptions are often intrinsically strong lines that can be detected at great heights 
above the stellar surface, giving us the ability to probe winds to much greater height. For 
the first time we could use the wings of Lya to measure hydrogen column densities directly 
for many lines of sight through the wind/chromosphere. The many lines of Fe II, likewise, 
gave measurements of excitation temperature (T exc ~ 5,000-12,500 K) and kinematics of the 
wind throughout much of the wind and upper chromosphere. 

Some of the results from IUE are conventional while others are surprising. Strengths 
and shapes of lines from different ionization stages of metals, most importantly iron, confirm 
the expectation that the metals are mostly singly-ionized throughout the chromospheres and 
winds of these stars. Likewise, the detection of the wings of Lya at height in chromospheres 
and rough agreement of the mass column densities derived with those from Fe II, means that 
H is primarily neutral throughout those parts of the wind we can sample. This contradicts 
the predictions of semi-empirical models in which H becomes completely ionized over the 
first several scale heights of the chromosphere, giving roughly a constant electron density 
in the line-emitting regions in spite of a marked decrease of total density (Judge 1990, p. 
290; § 3.2 below). This obsevation that H remains neutral thus places limits on permitted 
kinetic temperatures and local electron densities in the gas. Furthermore, wind models for 
31 Cyg give us some direct insight into the turbulence in chromospheres. Single-component 
models of the gas with no expansion effects, i.e., the sort of analyses used by Wilson & Abt, 
find Doppler widths of the order of 20 km s _1 , decidedly supersonic for gas with kinetic 
temperatures below 10,000 K. IUE observations show that at least some of this spread is 
caused by the differential expansion of the gas along the line of sight and need not be 
attributed to local turbulence. In fact, the IUE observations for 31 Cyg require a turbulence 
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< 15 km s _1 . 

Even in the classical optical analyses, the ionization of metals was lower than expected 
for a uniformly distributed gas and thereby implied clumping for the gas to achive enough 
electron density to maintain an observable population of trace neutral species, such as Fe I. 
UV observations confirm this result. Ionization throughout the wind is lower than expected 
from simplistic calculations of ionization equilibrium and implies clumping in the range 10- 
30x to achieve the inferred electron densities (n e ~ 0.2-1.5 x 10 9 cm -3 in the inner 150 R@ 
of the chromosphere). This is a complication well beyond most semi-empirical models of 
chromospheres. 

We can use an idealized description of the measured physical properties of the gas in the 
chromosphere and wind of 31 Cyg to test ideas about the structure of chromospheres and, 
ultimately, about wind mechanisms. These measurements are based on the most extensive 
set of observations ever obtained for a ( Aur binary (Eaton & Bell 1994). I have chosen to 
concentrate on 31 Cyg over the years because it has a much longer period than the other 
two classical ( Aur binaries (( Aur and 32 Cyg). This gives it a greater separation and 
less interaction between the wind and B star, a big advantage in many analyses. The wind- 
scattered emission lines seen in total eclipse, for instance, are weaker than in the other two 
stars (Eaton 1992). I am adopting the following properties for the 31 Cyg system: D = 
473 pc; R K = 197 R Q , M K = 11.7 M , M B = 7.1 M , Voo = 90 km s" 1 , ilS" = 3.0xl0- 8 
Mq yr _1 (Eaton 1993c; Eaton & Bell 1994). We should note that, in observing the more 
complicated ( Aur, Baade et al. (1996) probably had a better case than 31 Cyg for applying 
their wind-scattering analysis. Furthermore, inasmuch as the models for 31 Cyg and ( Aur 
give essentially the same results in terms of velocity structure/extension of the wind and 
wind temperatures, we can be confident that the results for either of them should be readily 
applicable to other stars. 

Table 2 gives the details of the empirical chromospheric/wind structure I am using 
to test wind models. Quantities listed are (1) radius (distance from center of the star in 
R & ), (2) expansion velocity, (3) excitation temperature, (4) a temperature to drive the wind 
thermally, (5) a clumping factor, CF, giving the inverse of the fraction of space actually filled 
with matter, (6) log(nn), the total hydrogen density, (7) log(n e ), the local electron density 
from an assumed constant ionization of H (3%) and amount of clumping, (8) log(n c ) for the 
variable ionization developed in § 3.2 and with the assumed clumping, but limited to 10% 
in the outer chromosphere, and (9) f equ , a velocity derived from equipartition between mass 
motion and internal thermal energy of the gas (see § 5.2). 

Although the electron densities measured from photoionization balance imply significant 
clumping of the gas, they are rather crude and may be systematically wrong. An independent 
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way to estimate this degree of clumping is to look at the difference between the pressureqj 
required to drive the observed wind and the gas pressures that would be available in the 
wind if it were not clumped. This approach works because the gas pressures (~ pT gas //i) 
would be in equilibrium with whatever pressure drives the wind. We do this by determining 
a temperature structure to drive the observed wind thermally, as though these stars had a 
coronal, i.e., generalized Parker- type, wind (see Lamers & Cassinelli 1999, Chapters 4&5), 
then comparing those temperatures to the observed temperatures, point by point, assuming 
the excitation temperatures we have measured approximate the electron temperature. This 
last assumption seems reasonable because the excitation temperatures are similar to temper- 
atures derived for semi-empirical chromospheric models of single stars, at least in the inner 
parts of the wind, but we will test it in later sections. To reiterate, if the gas is clumped, 
thereby bifurcating into dense clumps and an unspecified interclump medium, there must be 
a pressure in the interclump medium maintaining the clumping. We will assume that this 
is the pressure driving the wind and that it somehow breaks the wind up into microscopic 
clumps and maintains them. This is fundamentally different from the approach of traditional 
wave models in which the gas is uniform and the waves act on it continuously. 

The temperature structure to produce the pressure gradient required to drive the wind 
of 31 Cyg thermally, which we have derived from a thermal-wind model, is given roughly by 
the equations 

T th erm = 20, 000K + 75, 000K(r - R*)/R* for (r - < 15Oi? , (2a) 

and 

Ttherm = 95, 000K - 75, 000K(r - 3R*)/(15R*) for (r - R # ) > 2R # , (2b) 

where T therm is thermal (electron) temperature for a uniformly distributed gas, r is the radius, 
and .R*=197 R & is the photospheric radius of the star. I don't claim that this structure is 
a rigorous determination of the thermal-pressure profile, but that it is adequate for giving 
a good idea of the pressures required throughout the wind. Figure 1 shows the velocity 
structure calculated for this adopted temperature structure of Eq. 2; Figure 2, the amount 
of clumping implied. 

A way to test this sort of model is to calculate the emission expected from it. I have 
done this in two ways, first, with a spherically symmetrical model that calculates the line 



1 Strictly speaking, the wind is driven through the pressure gradient, which appears in the hydrodynamical 
momentum equation (F = ma), but a given pressure distribution, the higher level abstraction I am specifying, 
implies a pressure gradient, even in the simplistic case of an isothermal wind. 
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strengths, line profiles, and chromospheric mass column densities specifically for 31 Cyg. 
The second way uses a traditional plane-parallel model for ( Aur I developed in 1992-1995 
with PANDORA (Vernazza et al. 1973, Avrett & Loeser 1992) to explore the ionization 
structure and effect of clumping in these stars. 

3.1. Tests of a Spherically Symmetrical Model for 31 Cyg 

The spherical model for 31 Cyg has the distribution of gas given in Table 2. The 
velocity structure and distribution of mass density (given as nH=p/1.4mn) come from fitting 
spectra of shell lines in atmospheric-eclipse spectra; they should be fairly reliable. Excitation 
temperature comes from fitting the excitation of Fe II in these atmospheric-eclipse spectra; 
they should be reliably measured but subject to uncertainty about their meaning. Electron 
densities depend on level of ionization of H and degree of clumping. We can make educated 
speculations about these properties as follows. Since H seems to be neutral observationally, 
we might expect the ionization to be < 10%. We will assume it constant with height, 
following the thoughts of Judge (1990, § Illb), and of the order of 1-5%, and adopt a value 
of 3% for the sake of a first-order model. In that case, the gas must be clumped even to 
approach the canonical C II] electron density almost anywhere in the chromosphere. 

3.1.1. Line Emission 

For optically-thin lines, we can calculate the emission with the standard physics given 
by Osterbrock (1974), for example, with assumptions about ionization and chemical abun- 
dances. In this approximation, emissivity [ergs cm _3 s~ 1 ] is just 

e„ = 8.63 x 10" 6 fi lj2 n i n e exp(-x/A;T)T^ a5 cjr 1 x (3) 

where n^ is the density of emitting ions, Q is the collision strength, x is the excitation energy 
of the upper level, and lo\ is the statistical weight of the lower level. In this approximation, 
emission is proportional to collision strength, and we have incorporated collision strengths 
for C II] A2325 (Q = 0.830 and 1.66, for transitions out of the two ground-state levels, with 
half the emission going into the dominant 2325. 4A line) from Blum & Pradhan (1992) and 
for Al II] A2669 (fi = 3.3) from Tayal et al. (1984). Abundances of C and Al were A C /A H 
= 2.0 x 10~ 4 and A M /A U = 3.0 x 10~ 6 (Eaton & Bell 1994, § 2.1). To determine the total 
emission in an optically thin line as a function of radius on the sky, we sum the emissivity 
along a ray through the atmosphere, with integration limits determined by whether or not 
it intersects the star. To get a profile for the line, we define a velocity scale and map the 
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emissivity profile of the gas at each point in the atmosphere onto it, with allowance for the 
atmospheric expansion, t> exp , an arbitrary (systematic) velocity along the line from the center 
of the star, v syst , and a spectrum of Gaussian turbulence, v . 

The first thing we note from these models is the effect of the extended spherical at- 
mosphere on the emission strength and line profiles. Optically thin lines formed in such a 
structure will be highly limb brightened and will combine contributions from the face of the 
star and from a halo beyond its limb in similar proportions. Figure 3a illustrates this effect 
for our 31 Cyg model for an emitter assumed to exist in the same ionization stage throughout 
the wind (n c as in Table 2, Column 7). The emission peaks somewhat more to the limb than 
mass because electron density drops rapidly with height in these atmospheres, even with 
the assumed clumping. An emitter that is depleted in the inner chromosphere, as C II] may 
be, would peak even further beyond the limb. This is shown in the figure as a dashed line 
calculated for emission only at temperatures in the model above 5200 K. Even so, most of 
the emission would come from within the first < 0.15 i?* of the chromosphere/wind. Strong 
resonance lines, on the other hand, should be limb darkened. Scattering in an extended 
atmosphere would give a highly forward-peaked source function (e.g., Cassinelli & Hummer 
1971), the reason for the core-halo profile of the discs of WR stars, and, for resonance lines in 
chromospheres, the large optical depths mean the escape probability for photons migrating 
through the damping wings would be much smaller for the radial direction than for other 
lines of sight. 

Simple calculations like those in Figure 3a cannot represent the conditions in actual 
chromospheres because they give fluxes much higher than observed. Table 1 gives fluxes 
calculated for three cases, (1) uniform 3% ionization of H, (2) variable ionization of H, 
limited to 10%, and (3) variable ionization of H, limited to 3%. The model with uniform H 
ionization arbitrarily set at 3% gives an Al II] flux too high by a factor of ten and a C II] 
flux at least as bad. The model with variable ionization (Table 2, Column 8), on the other 
hand, gives Al II] flux high by ~50% and C II] high by slightly more, but the high electron 
density in the outer winds of these models, combined with single ionization of metals, gives 
much larger line broadening from differential expansion of the atmosphere. These models 
also have problems with electron density. If we average n e over emission of C II], the values 
are either too high, as for the calculation for the unrealistic uniform 3% ionization of H, or 
somewhat low, as for the realistic ionization of H limited to 3% or even 10%. Figure 3b 
shows the intensity profiles for the model with variable ionization. The effect of the 10% 
ionization in the wind is obvious, with much of the emission in the wind coming from the 
hot ionized outer parts. These parts, however, have electron densities much too low (7.7 in 
the log) to give the C II] line ratios, at least the ones seen in single stars, and they produce 
C II] profiles much broader than observed. They are clearly inappropriate for single stars 
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unless C becomes doubly ionized in the outer wind or unless the temperature in the outer 
wind is well below the measured excitation temperatures. 

3.1.2. Radio Continuum Emission 

A further test of these spherical models comes from the radio continuum observations 
of Drake & Linsky (1986) and Drake et al. (1987), who detected a number of bright cool 
supergiants at the 0.1 mJy level (or 10~ 27 ergs cm~ 2 s _1 Hz _1 ). Their measurement for 31 Cyg 
(0.36 ± 0.07 mJy) gives a check on the consistency of the electron densities and temperatures 
we are assuming. The basic idea (see Harper et al. 2005) is to integrate the radio source 
function, S v = B u , along rays through the wind and then sum the resulting intensities over 
the stellar disc to get a kind of luminosity [ergs s _1 ster _1 Hz -1 ]. This quantity would then 
be converted to flux at the Earth by multiplying it by the solid angle of 1 cm 2 at the star's 
distance, D. The dominant free-free opacity from electron-proton interactions, corrected for 
stimulated emission, is just 

k v = ^.2Yln v n e T~ x ^v- 2AQ (4) 

(Harper et al. 2005, Eqn. 3), while the generally smaller contribution from neutral-hydrogen- 
electron interactions (H~) is roughly 

k v = 6.2 x 10- 35 n H n e (\/lcm) 2 T e (5) 

(Allen 1973, § 42 evaluated at 8000 K). For the values of T e ,CF,n e in Table 2 (Column 8 
for n c ), we calculate a flux at the Earth of 0.04 mJy at 6.2 cm, only about 10% of the 
amount measured by Drake et al. Figure 4 shows the calculated intensity profile, a slightly 
limb-brightened disc with a radius of ~ 3.8 R*. This failure to reproduce the flux seems to 
be a common problem in reconciling the radio emission of cool giants with their UV spectra. 
Carpenter et al. (1999), for instance, were unable to reproduce the 3.5-cm flux of A Vel with 
the rapid wind acceleration they found from UV lines. Harper et al. (2005), likewise, were 
unable to reproduce the radio spectrum of ( Aur with their model and found it hard even 
to get the flux level. 

To get emission as great as is observed for 31 Cyg, the outer wind must be essentially 
wholly ionized, as the large mass-loss rate of Drake et al. implies. If we let the outer wind 
become mostly ionized, we can at least approach the flux level observed. For instance, with 
H in the envelope 75% ionized above 8000 K, we get a flux of 0.22 mJy; for 100% ionization, 
0.28 mJy. Is this level of ionization reasonable? It seems consistent with observations of 
( Aur stars, especially 31 Cyg, but the evidence for single stars is ambiguous. Such high 
ionization should reveal itself through absorptions of doubly ionized species in spectra of 
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( Aur binaries and possibly through emission lines from these species in both the binaries 
and single supergiants. The broadening of these lines would be a measurement of the terminal 
velocity and turbulence of the outer wind. Let us look at the observations of these species 
in some actual stars. There is a weak emission feature at the position of Si III] A1892 in 
the eclipse spectra of 31 Cyg (IUE images SWP47335 and SWP47336). In addition, strong 
emission lines of higher ionization are seen in eclipse in all three classical ( Aur binaries. 
Multiplets Al III UV1 and Fe III UV34 in 31 Cyg have P-Cyg profiles (Bauer & Stencel 1989, 
Fig. 3) and are so strong they must be formed by scattering of light from the B star in the 
wind. The line widths at the bases of their profiles are roughly 170 km s -1 , about what one 
would expect for a wind at terminal velocity. Eaton & Bell found that these absorptions are 
formed primarily at velocities toward the B star and likely result from ionization of the outer 
wind by the B star, not by intrinsic radiation of the supergiant. Thus a large fraction of the 
outer wind is ionized in the 31 Cyg system. Since H requires less energy to ionize than the 
19 eV to ionize C + out of its metastable level, H is likely ionized, as well, in these regions of 
31 Cyg. 

Single stars may have lower ionization in their outer winds than the binaries with similar 
mass-loss rates, although there is evidence in P-Cyg profiles of wind lines that some of the 
metals are doubly ionized in many of them (see §4.2 below). Carpenter et al. (1999) detected 
lines of C III] and Si III] in A Vel with about the right broadening and strength to be formed 
in an ionized wind. However, they interpreted these lines as an indication of gas at ~ 
50,000 K, as though high ionization necessarily means high temperature. The problem 
with wind formation is that the computed line profiles are wrong. The calculated profiles 
are essentially square, having a central dip reflecting the absence of gas at low expansion 
velocities, while the observed profiles, although somewhat noisy, seem to have a central peak 
as though formed at least partly in turbulent gas at low expansion velocity. Therefore the 
emission lines of single stars unfortunately tell us nothing about the ionization of their outer 
winds. 



3.2. Calculations with PANDORA 

PANDORA is useful for exploring non-traditional models of chromospheres because 
it lets one specify an arbitrary distribution of turbulence to increase the scale height in 
hydrostatic equilibrium and to specify an arbitrary distribution of added electron density 
that can be used to simulate clumping. I have made exploratory calculations for three 
cases: (1) a chromosphere for e Gem (G8 lb) based on the model of Basri et al. (1981), (2) a 
generalized model of a ( Aur binary (i?*=150 R Q ) with temperatures and scale heights based 
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on my analyses of ( Aur (Eaton 1993a) and 32 Cyg (Eaton 1993b), and (3) models for a Tau 
to explore excitation of Ha in cool giants (Eaton 1995). In the paper about Ha, I explored 
the conditions necessary for exciting Ha and C II] in semi-empirical hydrostatic models, 
finding that a certain amount of clumping of the warmer gas was required to strengthen Ha 
and to have electron densities high enough to give reasonable C II] line ratios. 

These semi-empirical models can give us some insight into how the gas might be ionized 
at various depths in actual chromospheres. All the models, regardless of assumptions about 
scale height, clumping, and temperature, predict species such as Al will be singly ionized 
throughout the whole emitting atmosphere. For H and C, the calculated ionization increases 
with height. Figure 5 shows this effect for H. In it, I have parameterized height with the 
electron temperature, since, by assumption, temperature increases monotonically with height 
above a temperature minimum in all such semi-empirical models. The calculations show that 
log(n c /nn) increases from a minimum of —4.0 at ~ 2600 K, determined by ionization of the 
metals included in the calculation, to 0.0 (complete ionization of H) at ~ 9400 K. For C, 
complete ionization occurs near 5000 K in these models. The actual models have scale 
heights somewhat smaller than that of 31 Cyg, so the level of ionization could be higher in 
31 Cyg. 



3.3. A Further Question about the 31 Cyg Model 

One possible error in Eaton & Bell's analysis of 31 Cyg is the determination of the 
temperature from excitation of Fe II. We have assumed that the excitation temperatures 
measured are bona fide electron temperatures of the gas. This might be the case, in that 
most of these lines arise from excited metastable levels which would probably be in thermal 
equilibrium with the electrons, at least in the denser parts of the wind (e.g., Judge et al. 
1992). Harper et al. (2005) argued this point explicitly, while others deriving temperatures 
for ( Aur binaries (such as Eaton & Bell) have implicitly assumed it. However, the exci- 
tation temperatures in the outer parts of the winds may well be greater than the thermal 
temperatures, since these regions are bathed in the radiation of a B star, and since the elec- 
tron densities expected in these zones are lower than the critical density (~ 10 6 per Judge 
et al. 1992) for radiative processes to become important. In the model for 31 Cyg, these 
temperatures have very little effect on the emission at any wavelength, since the densities 
are very low. 

Contrariwise, all of these stars have rather high excitation temperatures in their outer 
chromospheres, regardless of the effective temperature of the B companion. We see this 
in Figure 6 which gives T cxc as a function of tangential mass column density through the 
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chromosphere. We see that the temperature rises to about 8500 K where the Balmer lines 
become optically thin. This agrees roughly with theoretical calculations for supergiants. 

4. IMPLICATIONS FOR SINGLE STARS 

Here I shall attempt to fit the rich lore of space-dimensionless analyses of line emission 
and absorption of single stars into the context of our knowledge of the legitimately one- to 
two-dimensional knowledge of ( Aur binaries. We might expect the wind structure of the 
binaries to be essentially the same as that of single stars because the strengths of intrinsic 
emissions formed in the wind and chromosphere of ( Aur binaries are similar to those of single 
stars (Schroder et al. 1988; Eaton 1992; Harper et al. 2005, Fig. 6), and the Ha profiles, and 
their variation, seem to be the same as well (Eaton 1995; Eaton & Henry 1996). 

4.1. The Reality of Semi-Empirical Models 

Most detailed analyses of cool stars' chromospheres are based on semi-empirical models 
in which one posits a temperature rise through the outer atmosphere, calculates a density 
structure from hydrostatic equilibrium and radiative transfer, and then calculates the emer- 
gent spectrum from the physical conditions derived (e.g., Kelch et al. 1978; Basri et al. 1981; 
McMurray 1999). These models are inspired by the solar chromosphere, in which there is a 
temperature rise from a minimum, determined by the location where non-radiative heating 
begins to dominate heating by photospheric radiation, to a point where the chromosphere 
merges with a transition region of rapidly increasing temperature, heated by downward con- 
duction from a corona. The increase in temperature with height follows naturally from 
observations of the Sun and from the physical consideration that temperature should rise 
with decreasing density to keep emissivity from falling precipitately with height per Equation 
(1). Of course, there might be other ways to organize a chromosphere. For example, the 
material at different temperatures might be more intimately mixed throughout the chromo- 
sphere, as McMurry et al. (1999) contemplate to some degree. Harper et al. (2001) likewise 
have speculated about a truely radical reinterpretation of a Ori in which the chromospheric 
line emission comes entirely from hot inclusions in a generally cool neutral wind, although 
Harper et al. (2005) did not attempt to apply this radical approach to ( Aur. Indeed, tem- 
perature measurements for ( Aur binaries, at least to first order, confirm some sort of general 
temperature rise with height. Figure 6 shows the relation between excitation temperature 
and tangential mass column density through the atmospheres for four stars. The temper- 
ature rise is obvious. While this evidence does not preclude mixing some hot gas with the 
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bulk of the warm gas throughout the chromosphere, it shows that the bulk of the gas behaves 
sort of like the gas in these classical semi-empirical models. 

Another fundamental property of semi-empirical models is the density structure. In 
the calculated models, it usually comes from hydrostatic equilibrium between gravity and 
thermal motions of the emitting gas (T < 10,000 K). This is problematic in that there are 
obviously other sources of momentum flux in a typical chromosphere, such as the turbulence 
we see in profiles of emission lines, which will extend the atmosphere, changing its mass 
and electron density structure. Furthermore, the mere existence of turbulence implies some 
sort of clumping of the gas, which would necessarily change the local electron densities and 
ionization structure. Moreover, such effects would change the transfer of radiation and escape 
of photons from the chromosphere. 

A third property of many semi-empirical models is a precipitate temperature rise to 
coronal temperatures defining the top of the chromosphere. This is not a necessary feature, 
especially in the windy giants. Nevertheless, McMurry (1999) used such a rise for a Tau to fit 
the Lya profile and calculate emission from highly ionized species. However, the observations 
of C Aur binaries do not require such a rise and may not even allow it; the profiles of C III] 
and Si III] seem to require formation at least partly in the inner chromosphere/wind, and 
the arguments of McMurry et al. (1999, 2000) suggest there are other places to excite C IV. 

4.2. Terminal Velocities of the Winds 

One of the basic tenets of our understanding of the windy giants is the idea that their 
terminal velocities are much lower than the surface escape velocity (e.g., Hartmann & Mc- 
Gregor 1980; Judge 1992; Harper 1996). These terminal velocities are probably not as well 
known as we think they are, and they may well be much higher than generally thought, 
especially in the normal K giants like a Tau, as Judge (1992) argued in a provocative paper 
about wind energetics. This is hardly a new idea (see, e.g., Judge 1992; Ahmad & Stencel 
1988; Kuin & Ahmad 1989), but it is certainly worth discussing further in the context of 
supposed differences between single stars and binary components. For the classical ( Aur 
binaries, we base the terminal velocity on measurements of shell lines seen at all phases 
(e.g., Eaton 1993c). We should be suspicious of low terminal velocities since the metals in 
the outer parts of the winds might well be doubly ionized, as they are in calculations of 
Harper et al. (2005, Fig. 5) for ( Aur and § 3.2 above for single stars. Furthermore, the 
recombination time for H, given the electron densities in the outer parts of our model (~ 
10 3-5 cm -3 ), is of the order of 1-100 yr. Semi-empirical models for single stars also become 
fully ionized, at least in H, in these regions, if that actually means anything. So we really 
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don't know what these winds are doing in a vast volume of space before they recombine and 
form the shell lines. Occasionally more of the velocity structure may reveal itself through 
abnormally dense winds, as in AL Vel in 1992 (Eaton 1994, Fig. 8) and in A Vel in 1990 
(Mullan et al. 1998). Furthermore, we should realize that the winds must be sweeping up 
interstellar gas in these shells (cf. Wareing, Zijlstra, & O'Brien 2007). This means that the 
shell velocities are, if anything, likely to be lower than the true terminal velocities. 

Carpenter et al. (1999) admit that we probably do not see much of the velocity structure 
in the traditional shell lines in many stars, while they contend that they have seen it all in 
A Vel. However, it is not clear from the line profiles of Mg II, O I UV2, and Fe II UV1 that 
one sees it even in that case. The velocity structure for A Vel in 1994 seems to be ionization- 
bound, with the maximum wind velocity sampled limited by ionization of Fe + to Fe +2 . In an 
ionization-bound wind, the relatively high density behind the ionization front would allow 
somewhat weaker lines to become thick at their normal rate, while the strongest lines would 
quickly saturate at the velocity of the ionization front, little more than the velocity of those 
somewhat weaker lines, to give the sort of levelling off seen in the highest measured velocities. 
The edge velocities for all these strong lines in 1994 are around 50 km s _1 , much greater 
than the supposed terminal velocity of 33 km s _1 . Furthermore, in 1990 the star showed 
absorptions to at least 80 km s -1 (Mullan et al. 1998), which places its terminal velocity 
close to what we think 31 Cyg has, if the increase really did come from lower ionization in 
the outer wind. Also, in comparing a Tau and 7 Dra, Robinson et al. (1998) found terminal 
velocities of 30 and 67 km s _1 , respectively, for stars that otherwise seem to have similar 
atmospheric structure. All this evidence suggests terminal velocities several times as large as 
they often seem, with the apparent terminal velocity dependent on just how much material 
is being ejected into the flow at a given time. 

There is another problem in interpreting these shell lines (cf. edge velocities of A Vel) that 
would lead to errors in the terminal velocity. Carpenter et al. have interpreted the widths 
of Gaussians fit to the shell lines as a measure of turbulence. This could conceivably be the 
case for Fe II, but it is clearly inappropriate for at least some lines, for instance, for Mg II 
in 7 Cru (Carpenter et al. 1995; Fig. 10), where Mg II seems to have a significantly lower 
expansion velocity than even moderate Fe II lines. Such corrections for wind turbulence 
seem wrong, from both observational and theoretical considerations. Observationally, in 
models for scattering in shells (e.g., Baade et al. 1996), the turbulence drops with height. 
Theoretically, one would expect the turbulent energy to go into accelerating or heating the 
wind and be essentially damped out by the time the wind reaches its terminal velocity. 
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4.3. Turbulence of the Chromospheric Gas 

To investigate the nature of the turbulence, we may calculate the profiles of broadening 
for various velocity distributions for the model of 31 Cyg. The best fits seem to be for 
isotropic turbulence. In Figure 7, the dashed line shows the profile for isotropic turbulence 
giving roughly the line broadening of single supergiants. A Gaussian fits it quite well, as 
expected for the assumptions in the calculation. The solid curve, on the other hand, shows 
the effect of imposing a f syst =10 km s" 1 downward velocity on this same turbulence. The 
profile is still Gaussian at the level of the plot, but it is shifted to the red (3.8 km s _1 ) 
and broadened slightly more (AFWHM = 1.0 km s" 1 ) by the variable projection of the 
systematic velocity into the line of sight. 

Calculations for anisotropic turbulence all give non-Gaussian profiles to some extent. 
Some of these are shown in Figure 8. Radial turbulence is especially bad in this regard 
(Figs. 8a, c) because the strongest emitting regions are near the stellar limb, where the 
turbulent motions would be mostly perpendicular to the line of sight. Tangential turbulence 
is considerably better but still decidedly non-Gaussian (Figs. 8b, d), with enhanced wings not 
seen in profiles of actual stars (see, e.g., Carpenter et al. 1991, Fig. 3b). Figures 8c and 8d, 
however, show that elliptical anisotropic turbulence, in which one component is much larger 
than the other, is considerably closer to the observed Gaussian shape, and the departures 
from the observed Gaussian profiles would be somewhat reduced by convolving the calculated 
profiles with a Gaussian point-spread function for a spectrograph, which I have not done. 
However, these profiles would still have overly broad wings not seen in actual stars. I chose 
the nature of elliptical turbulence to reflect an isotropic turbulence of 5 km s -1 from perhaps 
thermal motion and 20 km s _1 of either radial or tangential macro/micro-turbulent motions. 

Carpenter & Robinson (1997) and Robinson et al. (1998) have cleverly fit the broadened 
wings of C II] lines with an anisotropic turbulence somehow averaged over a uniform disc, 
but that model is clearly inappropriate in light of the large limb effects predicted by ( Aur 
binaries (Figure 3). In fact, a calculation with our model, but limited to gas within the limb 
of the star, gives the awful non-Gaussian profiles seen in Figure 8e. 

Given the distribution of mass in the model and the resulting limb brightening, the 
Gaussian profiles of intrinsic lines imply a fairly isotropic turbulence. They certainly do not 
allow mostly radial turbulence, and they make it unlikely that the turbulence is perponder- 
antly tangential. 
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4.4. Electron Density and Interpretation of the C II] Multiplet 

We have hardly any actual knowledge of the electron densities in the winds/ chromospheres 
of any of these stars. What we do know is how dense the regions emitting the bulk of C II] 
are. Everything else must somehow come from a model calculation. 

One of the triumphs of applying physics to cool stars is the use of C II] to derive precise 
electron densities for a large number of objects. The kinematic analysis of Judge (1994), 
however, undermines the importance of this result. The redshift he detects in most stars 
is difficult to reconcile with other aspects of UV analyses, since it requires a ~ 10 km s _1 
downdraft of all the emitting gas or something more extreme if only part of the gas is 
participating. It would be convenient if this redshift were the result of a systematic error 
in the energy levels, and this may ultimately prove to be the cause of it. The energies 
of the upper levels of C II] are quoted by Moore (1970) to the nearest 0.1 cm -1 , about 1 
km s _1 , while the wavelengths are somehow known to one more significant figure. Likewise, 
the enhanced electron densities derived for the red side of these profiles seem to manifest 
themselves in only one of the line ratios, as though it might be affected by an unrecognized 
blend, but the dominant 2325. 4A line would have to be the one affected. 

However, if we accept the reality of the shift and enhancement, C II] must be telling us 
something profound about wind acceleration. C is two orders of magnitude more plentiful 
than Al, yet lines of the two are roughly the same strength in cool giants. Very little of the 
C, therefore, may be emitting C II], and it must be produced under rather special conditions 
to give the observed multiplet ratios of single stars. The broad wings of C II] can be fit by 
emission in the expanding gas of the upper chromosphere/lower wind with the broadening 
produced by differential projection of expansion into the line of sight. However, for that 
approach to work, the inner regions of the chromosphere must be neutral in C, even more 
so than implied by ionization balance in semi-empirical models such as those of § 3.2 or 
measured in 31 Cyg. Figure 8f shows such a calculation in which the gas below 5200 K 
emits no C II] . This calculation gives a reasonable flux, but it does not fit any of the other 
properties of C II]. The average electron density in the emitting gas (5xl0 7 ) is nowhere near 
the canonical 10 9 , even with the clumping. The width of the profile is somewhat more than 
expected for a real star (FWHM = 35 vs. 30 km s -1 for a weak line like Al II]), the peak is 
shifted blueward by —0.8 km s _1 , and the wings are probably too weak relative to the core. 
Clearly a more radical departure from the expected ionization/velocity structure is required. 

Looking at this more radically, we can get an idea of the conditions required in the 
31 Cyg model by extending our assumptions about ionization. Given the velocity struc- 
ture, observationally all the C would be inert below 7000-8000 K to weight the emission at 
significant expansion velocity enough to give the extended wings (giving the same sort of 
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effect as the anisotropic turbulence dismissed in § 4.3), with the bulk of emission coming 
from locally dense, probably rather cool regions at depth that may have peculiar velocities. 
Again, this means that C would be much less ionized in these atmospheres than predicted 
by semi-empirical models, with the C II] emission coming from only a moderately small frac- 
tion of the gas. However, the similarity of Doppler widths of C II] and other optically-thin 
lines such as Si II] and Al II] suggests that both are formed in gas with similar turbulence, 
perhaps in similar parts of the atmosphere. A possible source of the redshift of the line 
core is preferential excitation in the downward-facing parts of downward-moving elements in 
the turbulence distribution. This fraction of the gas would have the greatest relative speed 
with respect to the outward momentum source and might well be expected to be subjected 
to the most violent accelerations. Figures 9 show two experiments in simulating this ef- 
fect. In Figure 9a the emission comes from two sources, all the gas above 8500 K emitting 
with the velocity structure of Table 2 and 30 % of the gas at temperatures below 6000 K 
with a 10 km s" 1 downdraft. Figure 9b shows the emission for a model in which the gas 
above 8000 K is emitting with no downdraft and only the cooler gas with log(n c ) > 8.6 is 
emitting with a 10 km s" 1 downdraft. Both calculations give reasonable fluxes and electron 
densities averaged over the emitting gas (3.5 and 5.5 xlO 8 , respectively) high enough to 
begin to approximate the observed global values. Thus calculations with the most realistic 
chromospheric model we can muster are at least consistent with Judge's suggestions about 
C II]. 

This analysis has necessarily been speculative, although it is much better constrained 
by actual observations than any alternative analysis, e.g., one based on a traditional semi- 
empirical model, could have been. I think it is fair to conclude that (1) the broad wings of 
C II] likely come from formation in the lower wind (the broadening coming from differential 
expansion) and (2) C may not be ionized in the way expected from the standard non- 
LTE calculations of semi-empirical atmospheres but may be subject to some other source 
of ionization associated with the driving mechanism that picks out only part of the radial 
component of the turbulence distribution. In all these calculations the wings of the C II] line 
are broadened by formation in gas with significant expansion velocity projected into the line 
of sight. The asymmetry favoring the blue wing over the red results from having some of the 
redshifted gas blocked by the star. This effect is seen in actual stars, e.g., a Tau (Carpenter 
et al. 1991, Fig. 3) and 7 Cru (Carpenter et al. 1995, Fig. 2). Furthermore, as a caveat, 
we should note (1) that C would be depleated by the first dredgeup (e.g., Luck & Lambert 
1985), although this effect is presumably included in my adopted C abundance, and (2) that 
the observations of ( Aur binaries do not seem to show mostly neutral C at depth. One 
should also keep in mind that C II] may eventually destroy a basic assumption of this paper, 
and practically every other one written about chromospheres and winds of cool giants, that 
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the chromosphere and wind are one continuous structure spread evenly over the face of the 
star. 



4.5. More Highly Ionized Species 

The emissions of highly ionized species such as C III] and C IV do not fit readily into our 
model for 31 Cyg. The level of ionization required is incompatable with the temperatures 
and densities measured in the outer winds of ( Aur binaries as well as with the level of 
ionization actually measured in these winds. However, the clumping of the gas throughout the 
chromosphere/ wind gives us ample opportunity to incorporate more highly ionized material 
deep in the chromosphere and possible ways to excite it. We have assumed that the gas is 
confined to blobs that fill ~ 10% of space, surrounded essentially by a vacuum filled with 
an, as yet unspecified, source of pressure. The energy densities required to drive the wind 
through gradients of this pressure are certainly sufficiently high to allow higher ionization 
than the dominant levels we see in the shell lines. 



4.6. Velocity Structure from Fe II Self Reversals 

We may use our model for 31 Cyg to test the techniques used to deduce wind structure 
from shell lines in single stars. The problem here is that such analyses for single stars give 
winds accelerating much more rapidly than those in ( Aur binary components. I have used 
my 31 Cyg model to predict what the velocities of the shell components of strong Fe II lines 
in A Vel might be if it had the same velocity structure as 31 Cyg. This is the star with the 
best evidence for wind acceleration in Fe II self reversams, for which Carpenter et al. derived 
a wind with /3=0.9, t> 00 =29-33 km s _1 , and Iv[= 3xl0~ 9 M Q yr" 1 . In this simulation, I 
calculated the line-center radial optical depth into the atmosphere for the Fe II and Mg II 
lines measured by Carpenter et al. and determined the position and expansion velocity for 
r=2/3, as though the intrinsic line emission core comes from the center of the disc. Table 
3 gives the results, and Figure 10 displays them. These calculations assume the constraints 
of Table 2, with turbulence from Column 9. For the mass-loss rate of 31 Cyg, the lines 
become optically thick very high in the wind, absent ionization of Fe + to Fe +2 (Table 3, 
Column 5), but if we reduce the mass-loss rate somewhat to 9xl0~ 9 M Q yr _1 , to partially 
reflect Carpenter's lower value, and allow for ionization of Fe to Fe +2 for n e < 2xl0 6 (Table 
2, Column 8), we get velocities in much better agreement with observed values, as seen in 
Figure 10. Unfortunately we cannot apply this sort of analysis directly to the self-reversals 
of Fe II lines in 31 Cyg because the emision is formed not in the inner chromosphere, as in 
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single stars, but by scattering light from the B star in a huge shell surrounding the binary 
system. 

These results show that we can actually reproduce the sort of observations taken to 
indicate a rapid acceleration of the chromosphere in a single star with a model accelerating 
much more slowly. The secret for doing this is the realization that the observed absorption 
is limited by ionization of the outer wind as discussed in § 4.2. Therefore, there is really no 
compelling reason to believe the wind acceleration in the single windy giants is necessarily 
any faster than in binary components. 

There is evidence in occultations, as well as in the Ha profiles, that chromospheres/winds 
of single stars are much more extended than commonly thought. An occultation in Ha 
emission for the single M supergiant 119 Tau at M2 lb (White et al. 1982) and an image of 
a Ori at M2 lab (Hebden et al. 1987) both showed Ha emission, hence significant optical 
depths in the Lyman lines, out to at least twice the photospheric radii. Also, all of these 
cool supergiants show asymmetric Ha profiles indicating the line core is formed in the wind 
at least some of the time (Mallik 1993). Furthermore, in their model of the wind of a Ori, 
Harper et al. (2001) find a velocity /density structure giving even slower acceleration than 
we find in 31 Cyg or they found (Harper et al. 2005) in ( Aur. 

5. DISCUSSION 

The importance of ( Aur binaries is that they tell us just where the bulk of the gas is 
in the inner few stellar radii of the stars' chromospheres or winds. We can measure column 
densities directly and get at least a reasonable idea of how dense and hot it is throughout 
this critical zone where most of the energy is injected into a star's wind. Furthermore, with 
models like those of §3.2, we know that this gas must be emitting the intrinsic lines that we 
would easily see absent the B companions. The structure is therefore constrained in ways 
it simply cannot be in any analysis of a single star. I have subjected the model for 31 Cyg, 
and by extension models for other ( Aur binaries, to tests based on various sorts of emission 
from their chromospheres/winds. That such a model reproduces the intrinsic C II], Al II], and 
thermal radio emission to ~ 50% with a minimum of fiddling is an excellent test of the reality 
of the structure derived. Furthermore, the clumping inferred for such models provides an 
excellent opportunity (n.b., this really is fiddling) to accomodate tenuous material emitting 
lines of the more highly ionized species. This agreement is reassuring, inasmuch as most of 
the properties of the model are constrained by observations. The chief uncertainties involve 
ionization of H throughout the atmosphere, which is not measured well observationally, and 
ionization of H and the metals in the outer wind (T c > 9,000K) where many of the self 
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reversals of strong lines would probably be formed. 

5.1. Alfven-Wave Models 

There are several fairly extensive investigations of the conditions required for driving 
winds by Alfven waves that provide potentially testable predictions (Hartmann & McGregor 
1980; Hartmann & Avrett 1981; Holzen, Fla, & Leer 1983; Kuin & Ahmad 1989). For the 
most part these suffer from limitations of assuming the gas is effectively coupled to the fields 
(i.e., fully ionized), from not exploring the mechanisms of transfer of momentum from wave 
to gas realistically (although Holzer et al. did begin this process), and therefore from relying 
on somewhat nonphysical theories of how the waves are damped. 

There are a couple of ways of applying a wave model for driving a wind. In the more 
fundamental approach, one would use theories of how such waves are generated and how 
they interact with partially ionized media to predict the structure and properties of winds 
theoretically. Holzer et al. (1983, § VI) actually calculated some models of this sort. Unfor- 
tunately, the theories required to do that are difficult to apply, and a knowledge of magnetic 
fields of actual stars that would support such waves is lacking. The alternative is to form a 
semi-empirical model for wave-driven winds, somewhat like the chromospheric models I have 
been discussing, and determine what the properites it must have to fit the observed structure 
of a wind (viz., mass-loss rate, terminal velocity, velocity-density structure). A good example 
of this second approach is the analyses of ( Aur binaries by Kuin & Ahmad. Their models 
find that damping of the wave amplitude must decrease with height to fit observed velocity 
profiles, as one might well expect theoretically to keep the terminal velocity consistent with 
observation. Their models also give predictions of the level of turbulence in the chromo- 
sphere/wind by associating the lateral displacements of gas by such waves with the observed 
Doppler widths of gas in the atmospheres of these stars. Predicted Doppler widths, both by 
Kuin & Ahmad and Hartmann & Avrett seem to be larger than the observed turbulence, 
both in shell absorptions in ( Aur binaries and in the optically-thin emission lines of single 
windy giants. Furthermore, since Alfven waves are transverse, the "turbulence" would be 
anisotropic to first order. This prediction is at odds with the observation of isotropic tur- 
bulence, although it is based on the idealization of radial magnetic fields, while the actual 
topology might well be more complicated (see, e.g., Cassinelli et al. 1995). 

Kuin & Ahmad found from their semi-empirical models that the damping length for 
Alfven waves must increase with height. The most convincing mechanism for transferring 
energy from the wave to the gas, i.e., damping it, is so-called ion-neutral friction in which 
there is a phase lag between the wave's transferring momentum to ions and the ions' sub- 
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sequently transferring it to neutrals. Hartmann & McGregor discussed this mechanism, 
although they had no way of applying it a priori For waves with long periods, the transfer 
can be so rapid that the neutrals are effectively bound to the ions through elastic collisions, 
and there is little dissipation of wave energy or transfer of momentum. For high frequency, on 
the other hand, the neutrals cannot respond fast enough to the passage of a wave to partake 
in its displacements, and the ions just stir them up and dissipate wave motion as heat. Since 
the dissipated wave energy would go primarily into heating, this mechanism would be better 
for heating the chromosphere/ wind than driving its outflow. Most of the wave energy would 
be available for heating the wind but not for driving it, since the measured temperatures in 
these winds are much too low to drive them by thermal expansion (cf. § 3). Presumably the 
momentum of the wave would be transferred into wind motion. However, since the ratio of 
energy to momentum (E/p) goes as 1/v, with the Alfven speed generally larger than either a 
thermal or turbulent velocity, philosophically an attractive driving mechanism would trans- 
fer its energy to heat or some other mass motion as an intermediate stage. This is why it is 
so hard to drive winds with radiation pressure (E/p ~ 1/c). If all the waves had the same 
frequency, we would expect the momentum to be deposited in a narrow range of density, 
hence height, contrary to measurements of the acceleration of actual winds. 

The mechanism of ion-neutral friction would imply a damping length that decreased 
with height contrary to Kuin & Ahmad's result. However, it also implies that waves of 
different frequencies would be absorbed at different height and, therefore, that the spectrum 
of Alfven waves would determine the velocity profile of the wind. At this point our knowledge 
of the photospheric motions that might be exciting Alfven waves and the spectrum they 
would produce seems too sketchy to make any testable predictions about a wind's velocity 
structure. 



5.2. An Alternative Wind Mechanism 

Let us now take the liberty to speculate about a different way of driving the mass loss 
of these windy giants. We have developed here a picture of what conditions are required to 
drive the wind of one particularly well observed wind structure. Pressures required are an 
order of magnitude greater than those of the implied density/temperature structure. That 
the microturbulence required to fit line shapes and widths seems to be rather isotropic means 
that we likely are not simply seeing the effects of globally organized Alfven waves passing 
through the gas, as proposed by Hartmann & McGregor (1980) and implied by models of 
Airapetian et al. (2000), for example. There may be another way of supporting a wind with 
magentic fields, namely using chaotic fields emerging from the star and diffusing through 



-25 - 



the gas into space to drag the gas along with it and away from the star. Some form of this 
idea was implicit in our previous musings about the wind structure of 31 Cyg (Eaton & Bell 
1994, § 6), and Mullan et al. (1998) may have waved at it in passing. This very speculative 
picture is fundamentally different from the Alfven model in that the magnetic field being 
lost can constitute a moderate amount of luminosity. In the standard Alfven model, gas is 
driven ~ radially from the star along magnetic flux lines anchored permanently in the stellar 
surface. Here, magnetic flux would be lost at ~ the same rate as the gas and constitute a 
significant component of the energy loss in the wind through its adiabatic expansion. Such 
chaotic field would give a much more isotropic pressure, which would impress itself on the 
random velocities of the gas. We may estimate the effect by assuming the pressures driving 
the wind are in equipartition with the kinetic energy of random motions (turbulence) in the 
gas. For this condition the driving energy is ^NkT per degree of freedom (three of them), 
N being the number of particles in a random blob of mass M. The kinetic energy of the 
blob is |Mf^ qu (per degree of freedom), so that equipartition gives Mv^ = NkT, or 

v cqu = (k/<m >) a5 T t ° hc 5 rm = 0.08T t ° hc 5 rm (6) 

where we have taken the average mass per particle to be 1.3m^. This is to within a factor of 
^7 of the sound velocity, but for our elevated artificial temperature. For values of the driving 
temperature, T t h C rm (Table 2, Column 4), we get the turbulent velocities given in Table 2, 
Column 9. These values are comparable to the line-of-sight random velocities measured in 
cool (super)giants, and this fact argues that the driving force must be able to produce the 
random motions observed. It is unlikely, therefore, to be global Alfven waves. Of course, 
the pressure of these turbulent motions is itself a major source of momentum and energy in 
extending the atmosphere and driving the wind. 

Energy input determines the magnetic fields required in this model, since the magnetic 
energy density must be greater than the energy per unit volume required to lift the mass out 
of the potential of the star, B 2 /87r > GM-^p/R*- For our model, this leads to a field strength 
of 25 Gauss. If we assume the energy loss is double the potential energy from kinetic energy 
of the wind and emission from it, the field strength increases only to 35 Gauss. 

This kind of driving has the advantage over Alfven waves of being able to begin to 
explain the variation of mass loss from star to star in a way related to stellar structure 
and evolution. In our simulation of driving with thermal profile, the expansion-velocity 
structure is determined by the temperature (i.e., energy-density) profile, while the mass-loss 
rate is arbitrary, determined by supplying enough energy at some M to maintain the energy- 
density profile. With magnetic-flux emergence as the driving mechanism, the magnetic 
energy corresponds to the infusion of thermal energy in the coronal model; mass loss rate, 
therefore, is proportional to the rate at which magnetic field emerges. This is very attractive 



-26- 



in that there are indications that the winds of the giants vary in response to changes that can 
re-excite dynamos in their cores. See Mullan & MacDonald (2003) for changes of mundane 
giants, and recall the shells episodically thrown off by pulsing AGB stars. Alfven waves, 
on the other hand, would most likely be excited by convection, therefore be proportional to 
luminosity, and be little affected by the strength of the passive magnetic fields serving as 
their medium. 

As an alternative to chaotic magnetic fields filling the voids in chromospheric gas, we 
may imagine Alfven waves trapped in the cavities between ionized blobs. These waves would 
have speeds approaching the speed of light as the density dropped, and would be reflected 
off the blobs if their frequencies were below the cyclotron frequency. The critical frequency 
would rise as the material became more highly compressed (denser). Eaton & Bell actually 
had this mechanism in mind as the driving pressure in such atmospheres. There must be a 
rich optics of these Alfven waves waiting to be discovered. 

6. SUMMARY 

I have constructed a model for the chromosphere and wind of 31 Cyg which is based on 
measurements of physical properties in the outer atmospheres of 31 Cug and other classical 
( Aur binaries. It goes beyond other such models in that it derives the poorly understood 
turbulence and clumping of the gas from the pressures driving the wind's expansion. It 
predicts emission of optically-thin lines and microwave continuum to within 50 % of observed 
values, excellent agreement in the circumstances. In this model, the momentum flux required 
to drive the wind determines the stratification of the chromosphere where intrinsic lines 
would be formed. That momentum flux gives much lower densities than the stratification 
from thermal momentum flux alone, with the consequence that the gas must be clumped 
to produce the observed flux in intrinsic lines. We also find that this model can reproduce 
most of the properties of single stars' chromospheric spectra and argue that the evidence for 
fundamental differences between single stars and these binary components is rather weak. 

We must keep in mind, however, that there are some inconsistencies in this picture. The 
ionization of C, for instance, is a problem. It cannot be mostly singly ionized, as it seems 
to be observationally, without giving fluxes much larger than observed, and the model for 
31 Cyg does not predict the redshifts seen in single supergiants. C II] multiplet ratios from 
the model likewise would not predict the large global electron densities found from C II] in 
single stars. However, there are also inconsistencies in the interpretation of C II] in the single 
stars themselves, as I discuss in § 4.4. 
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I dedicate this paper to John S. Mathis who seemed to have the visceral intuitive grasp 
of physics required to address the sort of problems I have discused in this paper. If I'd had 
the sense to take his advice to work in this area when I was in graduate school, I probably 
would have written it years ago. This research has been supported by NSF through grant 
HRD-9706268 and NASA through grant NCC5-511 to TSU. 
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Table 1. Line Strengths in Cool Giants and Models 



Star 


Spectral 


V 


f(Lya) 


f(2325) 


f(2669) 


EW(Ha) 


log(n e ) 




Type 


(B-V) 






hine/W 






a Boo 


Kl III 


-0.04 


>1.4xl0 3 


35.7 


83.1 


1.12 


9.7 






1.23 


>3.7xl0~ 2 


9.4xl0~ 4 


2.2xl0~ 3 






a Tau 


K4 III 


0.85 


>340 


31.3 


30.7 


1.12 


9.0 






1.44 


>2,0xl0~ 2 


1.9X10" 3 


1.88X10" 3 






f3 Gru 


M5 III 


2.13 


>50 


31.7 


28.8 




8.5 






1.57 


>9,7xl0~ 3 


6.2xl0~ 3 


5.6xl0~ 3 






A Vol 


K5 Ib-II 


2.21 






12 


1.52 


8.9 












2.5xl0~ 3 






31 Cyg 


K4 I 


3.79 




<5 


2.3 


1.50 








1.28 




<5xl0~ 3 


2.1 xl0~ 3 






32 Cyg 


K4 I 


3.98 




<5±100% 


3.4 


1.75 








1.52 




<5.3xl0~ 3 


3.6xl0~ 3 






C Aur 


K4 lb 


3.79 




<4±100% 


2.2 


1.57 








1.22 




<3.7xl0~ 3 


2.0xl0~ 3 






Model 1 








56 


28 




9.8 


(3% ioniz) 
















Model 2 








9.1 


3.6 




8.5 


(var <10%) 
















Model 3 








6.7 


3.0 




8.6 


(var <3%) 

















Note. — Line fluxes at the Earth are in 10 13 ergs cm 2 s 1 . f(line)/fy is in A and assumes iy = 
3.65x10 ergs cm" - 2 s" 1 A- 1 at V = 0.0. Values of log(n e ) in Col. 8 come from C II] multiplet ratios as 
noted in the text. 
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Table 2. 


Details of Model for 31 Cyg 
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(3) 
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log(n H ) 
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log(n e ) 
(3%) 
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log(ra e ) 
(var}» 
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l>cqu 
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4.13 


4.66 
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3.79 


5.49 


4.55 
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12500. 
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5.65 
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19.0 
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ionization as in Figure 5 but limtcd to 10% in the outer chromosphere. 
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Table 3. Calculated Expansion Velocties for Shell Lines 



Multiplet 


A 


log(LSF) 


,.obs 
v XVel 


q ,calc 
U 31Cyg 


n.calc 
"Wei 




(A) 


©6000K 


(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 


Fe II UV1 


2599.40 


0.86 


-35.0 


-83.4 


-38.0 


Mg II UV1 


2795.52 


0.66 


-31.6 


-83.8 


-38.0 


Mg II UV1 


2802.70 


0.36 


-33.1 


-83.7 


-38.0 


Fe II UV1 


2598.37 


0.35 


-33.4 


-82.4 


-34.7 


Fe II UV1 


2585.88 


0.30 


-34.7 


-82.2 


-33.9 


Fe II UV1 


2607.09 


0.25 


-34.1 


-82.0 


-33.9 


Fe II UV3 


2332.80 


0.19 


-32.2 


-81.8 


-33.0 


Fe II UV62 


2755.73 


0.04 


-32.0 


-81.0 


-34.7 


Fe II UV3 


2364.83 


-0.01 


-30.0 


-80.6 


-30.5 


Fe II UV1 


2625.66 


-0.05 


-33.7 


-80.3 


-30.1 


Fe II UV32 


2739.55 


-0.05 


-27.2 


-80.3 


-33.9 


Fe II UV3 


2338.01 


-0.06 


-32.0 


-80.3 


-30.5 


Fe II UV1 


2617.62 


-0.16 


-33.4 


-79.4 


-28.4 


Fe II UV35 


2362.02 


-0.62 


-17.8 


-69.8 


-20.3 


Fe II UV35 


2331.30 


-0.75 


-29.0 


-64.0 


-16.6 


Fe II UV35 
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Fe II UV64 
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2591.52 


-0.93 
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Fig. 1. — Coronal wind model to fit the derived velocity structure of 31 Cyg. The dotted 
line is the velocity profile derived for 31 Cyg by Eaton & Bell (1994), while the solid curve 
is the velocity profile calculated with the temperature profile of Equations 2. Other curves 
are the local escape velocity (dashed) and sound speed for the velocities in Equations 2 
(dot-dashed). The circle shows the sonic point in this model. 
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Fig. 2. — Inferred clumping for the 31 Cyg model. Solid curve is the calculated clumping 
factor, CF, and the dashed line is the equipartition velocity defined by Equation 7. The 
atmospheric expansion velocity, v cxp is shown for comparison. 




-37- 





Optically-Thin Lines 


i i i l i i i i 
Cumulative 

■1 - ' 
■' 1 
\ 
l| 
'l 


3% ioniz of H 




h 

\ 

\ Emission@T>5200K 




\ 

\ '-.Moss _ 


| Emission 





100 200 300 

Radius on Sky, 5 (R Sun ) 



III 

Optically-Thin Lines 


Variable Ioniz. 


1 
ii 
n 
h 
1 1 

(/' \ Cumulative 
l N ■• 




Emission@T>5200K 




\ Emission 





100 200 300 

Radius on Sky, 5 (R Sun ) 



Fig. 3. — Calculated intensity of emission over the stellar disc for optically-thin emission 
lines, (left) This figure assumes the emitter is in the same ionization stage throughout the 
chromosphere, n e in Column 7 of Table2. (right) This figure is for the more realistic level of 
ionization (n e in Column 8 of Table2). 




Fig. 4. — Calculated intensity of 6.17-cm radio emission over the stellar disc for two models: 
dashed curve for the model with variable ionization limited to 10% (n c in Column 8 of Table 
2) and solid curve for the same model with 75% ionization of H in the outer atmosphere. 
This latter distribution would project a disc of 38 mas and have a spectral index a = 1.08, 
comparable to the few values measured by Drake & Linsky (1986). 
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Fig. 5. — Calculated ionization in three semi-empirical models for cool giant stars. The 
graph gives the ratio of electron density to total H density as a function of temperature. At 
low temperature, where ionization of H is low, the electron density is dominated by metals. 
At high temperature H becomes completely ionized, in these models, if not in actual stars. 
The region in which most of the chromospheric lines are emitted has a H ionization of a 
few percent. The solid curve shows the ralation I have adopted to relate electron density to 
excitation temperature in the calculations with variable ionization of H. Discrete symbols 
represent calculations for three models: dots, a Tau (Eaton 1995: Table 4, Model t8), circles, 
e Gem, and asterisks, ( Aur. 
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Fig. 6. — Excitation temperature, T cxc , vs. tangential mass column density, J pdx, as mea- 
sured in four ( Aur systems. Plusses are for ( Aur, asterisks for 32 Cyg, circles for 31 Cyg, 
and dots for 22 Vul. 
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Fig. 7. — Line profiles for the spherical model for two typical ions. The dashed curve shows 
the calculated profile for a line like Al II] excited throughout the chromosphere with an 
isotropic turbulence of 18 km s _1 superimposed on the expansion of the chromosphere/wind. 
This profile is fit to within the resolution of the plot by a Gaussian with FWHM = 30.3 
km s^ 1 . The solid curve shows the effect of superimposing a global 10 km s _1 downdraft on 
this profile. It is again fit with a Gaussian (no C II]-like wings) but shifted 3.8 km s _1 to 
the red and broadened slightly to FWHM = 31.3 km s _1 . 
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Fig. 8. The effects of various types of anisotropic turbulence on calculated line profiles. In these figures, the solid curve 

is the calculated profile, and the dashed curve is a Gaussian fit. For the more extreme profiles to the left, this Gaussian is fit to 
the whole profile; for the subtler profiles to the right, the Gaussian is fit to inner part (Intensity>0.5). Panel a) shows the really 
awful effect of purely radial turbulence, while Panel b) shows the more subtle effect of purely tangential turbulence. Panels 
c) and d) show, respectively, the effect of 20 km s _1 of radial or tangential turbulence combined with 5 km s _1 of isotropic 
turbulence. Panel e) is an experiment with confining radial and tangential turbulence to only the disc. The radial distribution 
is noticably different than Panel a) because most of the emission in Panel a) comes from beyond the edge of the disc. Panel f) 
shows the effect of isotropic turbulence with only the gas hotter than 5200 K emitting. 
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Fig. 9. — An attempt to fit the peculiar shapes and shifts of C II] emission in windy giants. 
Calculated profiles are the solid curves; Gaussians fit to the inner parts of these profiles 
(I>0.5) are dashed curves. In both cases gas in the deeper parts of the chromosphere (cooler 
than some threshold value) was given a 10 km s _1 downward velocity. Panel a) at left shows 
the effect of suppressing all the emission at temperatures between 6000 K and 8500 K while 
letting 30% of the gas below 6000 K emit. Panel b) at right shows the effect suppressing all 
emission from gas cooler than 8000 K except for gas with electron densities above log(n e ) 
> 8.6. In both cases extra broadening in the wings comes from emission from gas with 
significant expansion velocity projected into the line of sight. 
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Fig. 10. — A comparison of shell velocities for A Vel with calculations for the 31 Cyg model 
under assumptions about mass loss and ionization. Measured values come from Carpenter et 
al. (1999, Table 3). The circles are two lines for which the measured values seemed discrepant 
in Carpenter's paper, one even falling beyond the limits of his Fig. 9. Calculated velocities 
assume i\5~= 9xl0~ 9 M© yr _1 and that Fe is doubly ionized at n c < 2xl0 6 . 



